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We construct charmonium-like and bottomonium-like molecular interpolating currents with quan¬ 
tum numbers J PC = l + ~ in a systematic way, including both color singlet-singlet and color octet- 
octet structures. Using these interpolating currents, we calculate two-point correlation functions and 
perform QCD sum rule analyses to obtain mass spectra of the charmonium-like and bottomonium- 
like molecular states. Masses of the charmonium-like qccq molecular states for these various currents 
are extracted in the range 3.85-4.22 GeV, which are in good agreement with observed masses of 
the Z c resonances. Our numerical results suggest a possible landscape of hadronic molecule inter¬ 
pretations of the newly-observed Z c states. Mass spectra of the bottomonium-like qbbq molecular 
states are similarly obtained in the range 9.92-10.48 GeV, which support the interpretation of the 
(10610) meson as a molecular state within theoretical uncertainties. Possible decay channels of 
these molecular states are also discussed. 
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I. INTRODUCTION 

To date, there are eight members in the family of the electrically charged states: Z c ( 3900), Z c ( 4020), Z- v (4050), 
Z 2 ( 4250b Z c ( 4200), Z( 4430) and Zf,(10610), Zb{ 10650) observed in decays into final states containing a pair of heavy 
quarks [lHll|. Being not conventional QQ states because of their charge, they must be exotic with minimal quark 
contents QQud. 

The first charged exotic state, Z( 4430) + was observed in the B meson decay process B° —> i/j(2S)tt + K~ by the Belle 
Collaboration [l| in 2007. Recently, the LHCb experiment repeated the Belle analysis and confirmed the existence 
of Z( 4430) + with J p = 1 + [§]■ The broad doubly peaked structure Zi(4050) + and i?2(4250) + are resonances in the 
TT+Xci channel, which were found by the Belle Collaboration |3j in 2008. In 2013, the BESIII Collaboration reported 
Z c (3900) + in the process of V(4260) J/tfnr + TT~ 0, which was confirmed later by Belle 0 and CLEO data [f|. 
The BESIII Collaboration also observed the Z c (4020) resonance in the e + e - —> TT +i n~h c and e + e _ —> (D*D*)^^ 
processes 0,Q- Very recently, two new charged charmonium-like resonances Z c (4200) + [f>j and Z c ( 4050) [l0| were 
observed by the Belle Collaboration in the processes of B° —> Jif'K~K + and e + e _ —> n + n~i/j(2S), respectively. 
For the charged bottomonium-like states, Zb(10610) and Zb(10650) were reported bv the Belle Collaboration in the 
7r ± T(nS') and mass spectra in the T(5S) decay (llj| . One can consult Refs. fT^-ITil for recent reviews of these 
charged resonances. 

These exotic charged resonances are isovector states with quantum numbers I G J P = 1 + 1"*" while their neutral 
partners have charge-conjugation parity C = — 1. As four-quark states with quark contents ccud/bbud , these newly 
observed resonances were usually studied as hadron molecules and tetraquark states. These two hadron configurations 
are totally different. At the hadronic level, the hadron molecules are loosely bound states of two heavy mesons formed 
by the exchange of long-range color-singlet mesons. Tetraquarks are more compact four-body states which are generally 
bound by the QCD colored force between diquarks at the quark-gluon level. There are many theoretical studies on 
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these charged resonances; see Ref. 0 for a recent review. z? c (3900) + was interpreted as a DD* molecular state in 
Refs. 0-0. i? c (4020) + was speculated to be a D*D* molecular state in Refs. |18l - [20| . Z c ( 4200) + is much broader 
than other charged resonances in this family so that it was studied as a good candidate for a tetraquark state in 
Refs. 21, [23. It was also studied as a molecular state in Ref. [23j. Z(4430) + was described as a tetraquark state in 
Refs. 24-|26j and a D*D\ molecular state in Refs. [27l - l29j . Zb( 10610) and ^(10650) were interpreted as BB* and 
B*B* molecular states in Refs. Hi Hi. 

The mass spectra of the charmonium-likc and bottomonium-like tetraquark states were studied comprehensively 


in Refs. [3214351 |37j. In Ref. [35|, the masses of the charmonium-likc tetraquark states with I G J PC = 1 + 1 H were 
obtained from various currents in the range 4.0-4.2 GeV, which were consistent with the spectra of the charged Z c 
states. However, the molecular interpretations for these states are slightly more natural, especially for the Z c (3900) + 
and Z c ( 4020) + mesons which lie very close to the open-charm thresholds. In this work, we will study the mass spectra 
of the charmonium-like(qccg) and bottomonium-like(g66(j) molecular states with the quantum numbers J pc = l" 1 
using the approach of QCD sum rules Hsl PiTl . which is used to study the hadron properties of the lowest bound state. 
The masses of higher excited states are not easy to be calculate in QCD sum rules because their contributions are 
exponentially suppressed. However, there have been some attempts to study the orbitally excited nucleon (4ll - l43| . 
We try to explain the newly observed Z c and Zb states as molecular states and compare the difference between the 
molecular and tetraquark configurations. We note that we shall consider both the color singlet-singlet molecular 
structure and the color octet-octet “molecular” structure. 

There are many investigations of possible molecular states by using hadronic level Feynman diagrams, particularly 
investigations in the framework of the one boson exchange model 0- Generally speaking, this kind of study employs 
an effective Lagrangian to derive either the scattering amplitude or the effective potential. Besides the pion, quite 
a few other mesons are introduced which lead to many new coupling constants which have not been completely 
determined experimentally. Moreover, a form factor is always introduced at each vertex in order to suppress the 
high momentum exchange effect, which requires a new cutoff parameter. In other words, there exits some inherent 
uncertainties with the approach at the hadronic level. The QCD sum rule approach and the formalism at the hadronic 
level are complementary to each other. 

The paper is organized as follows. In Sect. II, we construct the molecular interpolating currents with J PC = l" 1 
for the Z c and Zb states. In Sect. Ill, we introduce the QCD sum rule formalism concisely and calculate the two-point 
correlation functions and spectral densities using these interpolating currents. We perform numerical analyses and 
extract the mass spectra and coupling constants for the charmonium-like and bottomonium-like molecular states in 
Sect. IV. In the last section, we summarize our results and discuss the possible decay modes for these molecular states. 


II. MOLECULAR INTERPOLATING CURRENTS FOR THE Z c AND Z b STATES 

There are two different types of four-quark operators: diquark-antidiquark type tetraquark fields and meson-meson 
type molecule fields. The former kind of operator is composed of a pair of diquark and antidiquark fields ([qq][qq]) 
while the latter one is composed of a pair of meson (or meson-like) fields (\qq \\qq]). The diquark-antidiquark type 
tetraquark fields have been constructed and studied systematically in Refs. [3J, |35|, [37| . Particularly, there are eight 
independent di qua rk-antidiquark tetraquark fields with J PC = D , which have been systematically constructed and 
studied in Ref. |35[. These eight tetraquark currents can be transformed into the combinations of other eight meson- 
meson type molecule currents by using Fierz transformations. In this paper, we systematically construct these eight 
meson-meson type molecule fields, and use them to study the charged Z p resonances as molecular states. 

The color structure of a molecule field([g<5][Q^]) can be written as 

(3 ® 3)[ 9 qj <8> (3 (g> 3 )[q^j = (1 ® 8)[ 9 q] ® (1 ® 8)[q^] 

= (1 ® 1) © (1 ® 8) © (8 ® 1) © (8 <g> 8) 

= 1 0 8 0 8 © (l © 8 © 8 0 10 0 TO © 27) . (1) 

The two color singlet structures in Eq. (JT)) come from the (l[ g gj <g> 1 [Qq]) and ® 8[Qg]) terms in the second step, 

respectively. In other words, the two mesonic fields [ qQ ] and [Qq] should have the same color structures to compose 
a color singlet molecular current. 

For (1 [qQ] <8 1[Q<?]) structure, we can obtain eight independent molecule fields with J p = 1 + by considering only 
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S-wave of the angular momentum between the two mesonic fields. Four of them are 

4 J, L = (Qa”f 5 Qa){QbJ/j.<lb), 

A l] L = (qaj a Qa){QbVa^qb), (2) 

Afi,L = 75 Qa){Qb&a.ii ( lb) •> 

while the other four can be obtained by performing the charge conjugation transform to these operators: 

= (da7 I^Q a){Qb'Y5Qb')i 

4 J,R = (9o7/j'75f3a)(Qb9b)) 

J 3 1]r = (qaCT afJ .'y5Qa)(Qbi a qb), (3) 

4fi,R = — {(}aO'ufj,Qa)(Qb'Y 75<?b)- 

In these expressions the subscripts a and b are color indices, and q and Q represent light quarks(u, d, s ) and heavy 
quarks (c, b ), respectively. 

Similarly, we can construct eight independent molecule fields belonging to (8 [ 9 qi ® 8[Qgj) color structure. Four of 
them are 

J^l = {qaisKbQbWci^KM, 

41]l = (qaKbQb)(QcWX n cd q d ), 

4*J,l = (qai a KbQb){Qca a ^Kdqdl ( 4 ) 

4%l = (qai a isKbQb)(Qc<j aiI .K d qd), 

while the other four can be similarly obtained by performing the charge conjugation transform to these operators: 

4J,R = ~(qal^abQb){Qcl5^c d qd)i 
4 %r = (qawKbQb)(QcXM, 

41]r = (qa^l 5 Xa b Qb)(Qc'y a Xc d q d ), (5) 

41 ]r = Aqa^KbQb){Qci a ibK d qd)- 

In these expressions A ra are eight Gell-Mann color matrices. 

These 16 molecule fields with J p = 1+ in Eqs. © © are independent, but they do not have definite charge- 
conjugation parities. We can use them to compose the molecular currents with definite charge-conjugation parities. 
The molecular currents with negative charge conjugation parity are 


J ( . 1} = J™ 


-J 


(i) 


r(8) _ r(8) _ n 
’’ J iu J iu.L 


8 ) 


'ifi u ifi,L °°i[i u ifi,L tJ ifi,R’> 
and the molecular currents with positive charge conjugation parity are 


T'f 1 ) _ n 1 -) i n ±J j'y°> _ no; i r 

J in.L ' J in..Ri J iu J in.L ' 


d 1 


rf 1 


'( 8 ) 


( 8 ) 


t( 8 ) 




ifi 




J ifi,R 


i= l,--- ,4, 

(6) 

,* = !,••• >4. 

(7) 


In this paper, we will study the Z c states by using the molecular currents constructed in Eq. © with quantum 
numbers J PC = l -1 

Jlfi = {qa'YbQa) {.Qb'y^qb) T {qa'Y/xQa) {Qb^bqb) i 

^2 fj, = {qaQa)(Qb'y^"/5qb) (da7/i75Qa)(Qb ( Zb)) 

J 3 ^ = {qa'Y Q a){Qb@aii'Ybqb) {qa&a/i'YbQ a) (Qb7 db); 

J4jj, (da7 75^a) (Qb^a^Qb') T (qa&otfiQa) (Qb'Y 'YbQb')i 

J i 8 J = {qalsKbQb){.Qc"1vK d q d ) + (ga7M A ab < 3b)(Qc75A” d g d ), (8) 

4 l ] = (qaK b Qb)(QcwK d q d ) - (Qa'r^KbQbKQAM 

J 3 S J = i.qal a X^ b Qb){Qc^ a ^bXZ d q d )-{qaCT a ^bXabQb){Qcl a X 7 i d q d ), 

4$ = (da7 a 75A^Qb)(Q c <7« M A^ gd ) + (g Q a QM A^g b )(Q c 7 a 75A^ g(i ), 
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in which j}*] — belong to the color structure (1 [ q Q] ® l[Qq]) while j[^ — jff belong to the color structure 
(8[ 9 g] ® 8[Qq])- The eight independent diquark-antidiquark tetraquark fields with J PC = l" 1 constructed in Ref. jUJ 
can be written as combinations of these eight meson-meson type molecule currents. Moreover, one can construct 
other eight “meson-meson” type molecule currents, having the color structure [q a Qb][Qb<la\- They can also be written 
as combinations of these eight meson-meson type molecule currents, having the color structures [q a Qa] [QbQb] and 
[qaX n ab Qb][QcX cd.qd\- I 11 general, there is no one to one correspondence between the current and the state. The 
independence of the currents means that if the physical state is a molecular state, it would be best to choose a 
molecular type of current so that it has a large overlap with the physical state. Similarly, it would be best to choose 
a tetraquark current for a tetraquark state. 

We note that the interpolating currents listed in Eq. ([5]) should contain the quark contents uccu + deed, to be neutral 
molecular currents of / = 1. Such molecular currents have quantum numbers I G J PG = l - * - 1 -1 and thus couple to 
neutral Z c states. The corresponding operators with the quark contents need or dccu can couple to charged Z c states. 
They altogether form isospin triplets. However, we will work in the SU(2) isospin symmetry without considering the 
effect of isospin breaking in this paper, i.e., we neglect instantons because we are in the vector channel, the masses of 
the up and down quarks and maintain isospin for the quark condensates (uu) = {dd) = {qq), ( uGu) = (dGd) = ( qGq). 
Accordingly, the QCD sum rules for any iso-triplet are the same. Moreover, the isoscalar molecular currents can also 
be obtained from Eq. © with the quark contents uccu — deed , and the sum rules for theses currents are also the same 
as those for the isospin triplet currents. Therefore, the same mass predictions would be obtained for the neutral and 
charged Z c states with I G J p ^ = l + l +( ^ and their isoscalar partner with I G J PG = 0 _ 1 H . This expectation is 

reasonable for these quarkonium-like states, for example, the neutral states Z®( 3900) ,(| and Z°( 4020) 136:] lie very 
close to their charged partner Z c (3900) + and Z c (4020) + , respectively. 


III. QCD SUM RULES FORMALISM 


With these currents constructed in Eq. ©, we can study the following two-point correlation function 

IV(g 2 ) = i J d i xe iq ' x {Q\Tyj^(x)Jl(Q)\ |0) 

= n(fl(^ + n'(, 2 ) (9) 

where n (q 2 ) and n'(g 2 ) are invariant functions related to spin-1 and spin-0 states, respectively. The two-point 
correlation function can be described at both the hadron and quark-gluon levels. At the hadron level, the correlation 
function has a dispersion relation representation 


n (<7 2 ) 


(Q 2 ) N 



Imn(s) 


s N (s — q 2 — ie) 


N -1 

ds + ^2 b r 

71=0 


(q 2 r 


( 10 ) 


in which b n are the unknown subtraction constants which can be removed by taking the Borel transform. The lower 
limit s< denotes a physical threshold. With this expression, one only needs to evaluate the imaginary part of the 
correlation function, which is much easier than the full calculation. The imaginary part of the correlation function 
is defined as the spectral function p(s) = ilmn(s), which is usually evaluated at the hadron level by inserting 
intermediate hadron states 

P{s) = J2S{s-ml){0\r]\n)(n\if\0} 

n 

= fxH s ~ m x ) + continuum , (11) 


where we have adopted the usual pole plus continuum parametrization in the second step. All the intermediate states 
\n) must have the same quantum numbers as the interpolating currents J^x). The lowest-lying resonance with hadron 
mass mx couples to the current J^{x) via 


(OIJ^I X)=f x e„, (12) 

in which fx is coupling constant and e M is the polarization vector (e • q = 0). 

At the quark-gluon level, we evaluate the correlation function and spectral density via the QCD operator product 
expansion (OPE) up to dimension-eight at the leading order of a s . The correlation function n(g 2 ) and spectral density 
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p(s) can be expressed in terms of quark and gluon fields. These results are compared with Eq. (1101) obtained at the 
hadron level to establish sum rules for hadron parameters, such as masses, magnetic moments and coupling constants 
of ground state hadrons. As mentioned above, we usually take Borel transform to the correlation functions at both 
the hadron level and quark-gluon level to remove the unknown constants in Eq. ca and suppress the continuum 
contributions. Using the spectral function defined in Eq. CD- the sum rules can be obtained as 

C k (s 0j M 2 b ) = [ ° dse- s / M *p(s)s k = f 2 x m 2 x k e- m x/ M B , (13) 

Js< 

where Mg is the Borel parameter and p(s) in the integral is the spectral density evaluated in QCD side. The upper 
integral limit Sq is the continuum threshold above which the contributions from the continuum and higher excited 
states can be approximated well by the spectral function. Finally, the hadron mass mx for the lowest-lying state is 
extracted as 


m x = 


I Ci ( sq , Ml) 

Co (so, Ml) 


(14) 


It is shown in Eqs. CD and CD that the extracted hadron mass mx is a function of the continuum threshold So 
and Borel mass Mg. One can perform the QCD sum rule analysis with these two equations. At the leading order in 
a s , the spectral density p(s) in Eq. (flp is evaluated up to dimension eight, including the perturbative term, quark 
condensate (qq), gluon condensate ( g 2 GG ), quark-gluon mixed condensate ( qg s a ■ Gq) 1 four-quark condensate ( qq) 2 
and the dimension eight condensate (qq)(qg s a ■ Gq) 


p{s) = f? ert (s) + p {qq) (s) + pi uu >{s) + + pW*>(s) + p 


(GG), 


■A Gq)i 


{<lGq) 


0). 


(15) 


( 8 ) 

To illustrate our numerical analysis, we use the current J as an example and show its spectral density in the 
following: 


p[ 8)pert (s) = 


1 


Pi 


( 8 )(qq) 


(s) = - 


192t r 6 

mQ(qq) 


roc-max r 

/ da 

•j Oim .in. J & 


Pmax (1 — a — /3)(1 + a + P) [(a + (3 )mg - aps\ 


dp 


a 3 P 3 

(1 — a — ft) [(a + P)mQ — aps\ [3 mg(a + (3) — 7a/3s] 
ap 2 

7toq(1 — a — /3) 2 (5 + a + p) [(a + P)m 2 g — a/3s] 
96a 2 /3 2 


(8 )(GG), s. 

Pi (s) = - 


67T 4 

(PsGG) 

288?r 6 


f 

J Otr 


da 


da 


dp 


dp< 


7toq( 1 + a + P) [(a + P)mg — a/3s] Wq( 1 — a — p) 2 \m 2 g(a + P) — 2a/3s] 

16a/3 a 3 

(1 — a — P) [rag(3 + a + /3) + 2a/3s] [(a + P)mg — a/3s] 1 


8a 2 p 


(8){qGq) 

Pi 


( S ) = 


mQ{qg s cr ■ Gq) 
144tt 4 


rf- 

/ da 

J OLm.iT,. J &r. 


dp 


(1 — a — p) [6mg(a + /3) — lla/3s] 17 [3mg(a + /3) — 5a/3s] 1 


2/3 


1 


„(8 )(qq) 2 
Pi 


J 8 )(qq)(qGq) 

Pi 


[s = 


(s) = 


4 m 2 Q (qq) 

9n 2 

(qq){qgs<7 ■ Gq) 
108tt 2 


1 — 4mg / s , 


da 


(16) 


\ 2Am Q : , r ^ , m 2 Q (5a 2 -6a + 3) \ 

J S ( s - m Q) + - Q(1 _ a)2 - 5 ( s “ m Q ) - aH ( s ~ m Q) } > 


a* 


in which frig = , S'(s — mg) = dS ( s ^ Q ) ; anc j H(s — frig) is a Heaviside step function. The integration limits 


are a n 


_ l-\/l-H/ 8 a _ 

— o ’ Hmin — 


= 1 — a, a„ 


= +y/ 2 and mQ is the heavy quark mass. 


2 5 f-'rrnn — ocs—m,Q ’ Mmax — ^ '-*5 max — 2 

We note that we have ignored the chirally suppressed terms with the light quark mass and adopted the factorization 
assumption of vacuum saturation for higher dimensional condensates (D = 6 and D = 8). The results for other 
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currents listed in Eq. ((Sj) are collected in Appendix. [A] From these expressions we can find that the D = 3 quark 
condensate ( qq ) and the D = 5 mixed condensate (qg s & • Gq) are both multiplied by the heavy quark mass rriQ , which 
are thus important power corrections to the correlation functions. 

IV. NUMERICAL ANALYSIS 

( 8 ) 

In this section we still use the current as an example and perform the numerical analysis. The following QCD 
parameters of quark masses and various condensates are used in our analysis (Til - Fill : 

m c {m c ) = (1.23 ±0.09) GeV, 
m b (m b ) = (4.20±0.07) GeV, 

(qq) = —(0.23 ± 0.03) 3 GeV 3 , 

(qg s a ■ Gq) = -M$(qq ), (17) 

Mg = (0.8 ± 0.2) GeV 2 , 

(g 2 s GG) = (0.88 ± 0.14) GeV 4 . 

Note that there is a minus sign implicitly included in the definition of the coupling constant g a . We use the running 
masses in the MS scheme for the charm and bottom quarks. 

As mentioned above, the extracted hadron mass mj in Eq. (1141) is a function of the continuum threshold sq and the 
Borel mass Mg, which are two vital parameters in QCD sum rule analyses. If the final result, nix, does not depend on 
these two free parameters, then the method of QCD sum rules would have perfect predictive power. However, reliable 
mass predictions are obtained when there is weak dependance on these parameters in a reasonable working regions. 
Principally, there are two criteria to find a Borel window(reasonable working region of Mb)'- the requirement of the 
OPE convergence results in a lower bound while the constraint of the pole contribution leads to an upper bound. 
At the same time, we will study the variation of the hadron mass mx with respect to the continuum threshold. An 
optimized value of the continuum threshold sq is chosen to minimize the dependence of the extracted hadron mass 
mx on the Borel mass Mg. 

In Eq. (1151) . the non-perturbative terms are evaluated up to dimension eight. After the numerical analysis, we find 
that the quark condensate (qq) and quark-gluon mixed condensate (qg s cr • Gq) are dominant power corrections while 
the contributions of other condensates are much smaller. Using the spectral density for the current in Eq. (fT71) . we 

show the contribution of the dimension eight condensate to the correlation function Il[ qq ^ qGq \Mg, so)/H1 u (Mg, sq) 
in Fig. [T] with so —> oo, in which the ratio decreases with respect to Mg. Accordingly, we require the dimension eight 
condensate contribution to be less than 5%, which results in a lower bound = 2.8 GeV 2 . 



M 2 b [GeV 2 ] 


(g\ 

FIG. 1: OPE convergence for the current 1 with so —► oo. 

To determine the upper bound on Mg, we define the pole contribution(PC) using the sum rules established in Eq. 

USD, 

Go (so, Mg) 

C 0 (OO, Mg) ’ 


PC(sq, Mg) 


( 18 ) 
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which represents the lowest-lying resonance contribution to the correlation function. The continuum threshold Sq 
is an important parameter to the pole contribution. We study the variation of mx with so in the left panel of 
Fig. El by varying the value of Mg from its lower bound. With different values of Mg, there curves intersect at 
so = 18 GeV 2 around which the variation of mx with Mg reaches it minimum. This is thus an optimized value of 
the continuum threshold so to study the pole contribution defined in Eq. m- Requiring that PC be larger than 
20%, we obtain an upper bound on the Borel mass M^ ax = 3.7GeV 2 . The Borel window is then determined to be 
2.8 GeV 2 < Mg < 3.7 GeV 2 with the threshold value so = 18 GeV 2 . 

In the right panel of Fig. El we show the variation of the hadron mass mx with respect to M|. The Borel 
window varies quickly for different value of so- One notes that the mass curves decrease significantly in the region 
Mg < M^ in while becoming quite stable inside the Borel windows. Finally, we can extract the hadron mass and the 
coupling constant 


to ( 8) = (3.90±0.12) GeV, (19) 

A 1 

/ (8) = (0.69 ± 0.21) x 10" 2 GeV 5 . (20) 

This value is consistent with the mass of Z c (3900), which implies the possible molecule interpretation of this new 
resonance. Here we would like to emphasize that in our calculations we have not used masses of heavy-light mesons 
as inputs, but simply note that the obtained value 3.9 GeV (as well as other mx listed in Table [T| is close to the 
threshold of two heavy-light mesons. Further studies are needed to understand whether there is an underlying reason 
for these results. 




ill; 


- 

-ill; 

■' " 


. Mj=2.8 GeV 2 



- - A/j=3.1 GeV 2 



— — Ml= 3.4 GeV 2 

■ - - - Mj=3J GeV 2 








3.5- 1 - 1 - 1 - 1 - 1 - 

10 13 16 19 22 25 28 

jo [GeV 2 ] 


4.6 

4.4 

? 4 - 2 
O 

g 4.0 
3.8 


s 0 =16GeV 2 
s 0 =18GeV 2 
s 0 =20 GeV 2 


- - I' 


3.6- 1 - 1 - 1 - 1 - 1 - 1 - 

1.8 2.2 2.6 3.0 3.4 3.8 4.2 4.6 

Mi [GeV 2 ] 


9 ( 8 ) 

FIG. 2: Variations of the charmonium-like molecule hadron mass mx with s o and M B for the current J) '. 

After performing similar numerical analyses for the other interpolating currents, we collect the extracted numerical 
results for the hadron masses and coupling constants in Tabic Q] The mass sum rules for the currents and 
are unstable and thus they do not give reliable mass predictions. For the current the lower bound of the Borel 
window is very small under the first (convergence) criterion. Although it leads to very good OPE convergence and 
broad Borel window, we need to consider the stability of the Borel curves, from which the lower bound of the Borel 
window is determined to be 3.3 GeV 2 . The situations for the currents and are very different, in which the 
lower bounds on Mg are bigger than their upper bounds, suggesting the OPE convergence is poor for them. By 
loosening the criterion of the OPE convergence and requiring the dimension eight contribution to be less than 10%, 
we can still obtain stable mass sum rules for jfj^f and and reliable mass predictions, as shown in Table Q] The 
error sources including the uncertainties of the various parameters in Eq. m and the continuum threshold sq are 
considered to obtain the errors for hadron masses and coupling constants. 



















Current 

so (GeV 2 ) 

Borel window (GeV 2 ) mx (GeV) 

fx (10- 2 GeV 5 ) 

jD) 

21 

3 . 3 - 

-4.5 

4.22 ±0.14 

0.57 ±0.16 

J 1 

J 2n 

- 



- 

- 

J 1 

20 

3 . 3 - 

-4.2 

4.04 ±0.12 

0.93 ±0.27 

t(!) 

j Al± 

20 

(3.0- 

-3.3) 

4.02 ±0.15 

0.35 ±0.13 

t(8) 

J li± 

18 

2.8- 

- 3.7 

3.90 ±0.12 

0.69 ±0.21 

7 ( 8 ) 

J 2fi 

- 



- 

- 

7 ( 8 ) 

J 3 [i 

18 

3.1 - 

-3.9 

3.85 ±0.11 

1.51 ±0.46 

7 ( 8 ) 

J 4n 

20 

(2.8- 

-3.1) 

4.03 ±0.18 

0.59 ±0.23 


TABLE I: Numerical results for the charmonium-like molecule states. 


In Table [TJ the extracted masses from the currents of color structure (8[ g gi (g> 8[Qgj) are about 3.85 — 4.03 GeV, 

which are slightly below the 4.02 — 4.22 GeV from the currents of color structure (l[ q Q] ® l[Qg])> although they 

( 8 ) ( 8 ) 

both lie precisely in the range of spectra of Z c states. The masses extracted from the currents and .7^/ are 
(3.90 ± 0.12) GeV and (3.85 ± 0.11) GeV respectively, which are clearly consistent with the mass of i? c (3900). The 
interpolating currents , jj^ and give hadron masses (4.04 ± 0.12) GeV, (4.02 ± 0.15) GeV and (4.03 ± 0.18) 
GeV respectively, which are in very close proximity to the masses of the Z c (4020) and Z c ( 4050) mesons, although the 
latter state is not confirmed to date. We note that these values are also in rough agreement with the mass of Z c (3900) 
state. However, one can find that it is better to chose the currents J\ and J V/ to fit the mass of Z c ( 3900) because 
these two currents have a larger overlap with the physical state. We can infer that Z c { 3900) has a structure well 
represented by the currents and ■ Last but not least, the current leads to a mass prediction (4.22 ±0.14) 


GeV, which is in good agreement with the mass of Z c ( 4200). The interpolating currents and J^ 1 

are constructed as DD *, D^Di, D*D* and D\Di molecular operators, respectively. Our results in Table [Q suggest a 
possible landscape of hadronic molecular interpretations of the charged and neutral Z c states. According to the above 
analysis, we suggest that the Z c ( 4200) meson to be a DD* state, while the Z c ( 4020), .Z c (4050) and .Z c (3900) mesons 
to be D*D* or D\Di states. The analysis of the current j! 2 ^ may imply that the stable D^Di molecular state does 
not exist. 


io.7 ptr 
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FIG. 3: Variations of the bottomonium-like molecule hadron mass mx with so and Mg for the current . 

Similarly, we can study bottomonium-like molecule states with J PC = l" 1 by taking mg = mb in the expression of 
the spectral density in Eq. 03 and AppendixjA] The bottomonium-like molecule system is similar to the charmonium- 
like system due to heavy quark symmetry. Under the same criteria, one finds that the Borel window for the qbbq 
system is much broader than the qccq system. This suggests a stricter limitation of the pole contribution, which is 
only required to be larger than 20% for the qccq systems. For the qbbq system with , we require the same OPE 
convergence criterion as the qccq system while modifying the requirement of the pole contribution to be larger than 
30%. The Borel window is obtained as 7.5 GeV 2 < Mg < 8.8 GeV 2 with the continuum threshold so = 108 GeV 2 . 
Using these values of the parameters, we show the variations of the hadron mass with respect to the Borel mass M|> 
and the threshold value sq in Fig. [3] The Borel curves are shown to be very stable and give reliable predictions of the 

























9 


hadron mass and coupling constant 


m v( s) = (9.93±0.15) GeV, 

A 1 b 

(21) 

/ (8) = (1.02 ± 0.30) x 10" 3 GeV 5 . 

b 

(22) 


After numerical analyses of all interpolating currents, we collect the numerical results for the qbbq states in Tabic HT1 
Similar to the charmonium-like system, there is no significant Borel window for the current under the above 
criteria. The Borel window (7.5 — 8.5) GeV 2 written in parenthesis is obtained by loosening the requirement of the 
OPE convergence to be 10%. However, the mass prediction under this Borel window is still reliable. The extracted 
mass for the current is about 10.48 ± 0.15 GeV, which is consistent with the mass of the 2f,(10610) meson within 
the error, supporting the BB* molecule interpretation for this state. 


Current 

s 0 (GeV z ) 

Borel window (GeV^) mx (GeV) 

fx (10- 3 GeV°) 

J 111 

121 

8 .0- 

- 11.0 

10.48 ±0.15 

1.45 ±0.19 

J (1) 

J 2 v 

- 


- 

- 

- 

A 1 ) 

J 3 /x 

113 

8.0 

- 9.4 

10.14 ± 0.15 

1.69 ±0.46 

jW 

117 

7.5 

- 8.3 

10.33 ±0.14 

0.58 ±0.34 

t(8) 

J 1 /x 

108 

7.5 

- 8.8 

9.93 ±0.15 

1.02 ±0.30 

t(8) 

J 2 n 

— 


— 

— 

— 

t(8) 

J 3 /x 

108 

7.8 

- 8.7 

9.92 ±0.15 

2.17 ±0.62 

t(8) 

119 

(7.5 

-8.5) 

10.46 ±0.14 

1.67 ±0.55 


TABLE II: Numerical results for the bottomonium-like molecule states. 


V. SUMMARY AND DISCUSSIONS 

To study the charged exotic Z c and Zb states, we constructed all the charmonium-like/bottomonium-like molecular 
interpolating currents with J PG = l+ _ , including both the singlet-singlet and octet-octet types of color structures. 
We calculated the two-point correlation functions and spectral densities for these qQQq operators. Within the 517(2) 
isospin symmetry, all the numerical results of hadron masses and coupling constants in Tables U and mi are suitable 
for the neutral and charged Z c states with I G J pi ' C ' ) = l + l + (~) and their isoscalar partner with I G J PC = 0 _ 1 H . 

At the leading order in a s , we calculated the two-point correlation functions and spectral densities up to dimension 
eight, including the perturbative term, the quark condensate (qq), the mixed condensate (qg s & ■ Gq), the gluon 
condensate (g 2 GG), the four-quark condensate (qq) 2 and the D = 8 condensate ( qq)(qg s <J ■ Gq). Being proportional to 
the heavy quark mass, the quark condensate (qq) is the dominant power correction to the correlation function while 
the mixed condensate (qg s & ■ Gq) also gives an important contribution. After performing the numerical analyses, we 
obtain reliable mass predictions 4.02 — 4.22 GeV for the color singlet-singlet charmonium-like molecular states and 
3.85 — 4.03 GeV for the color octet-octet ones. This mass spectrum of the qccq states is precisely consistent with the 
masses of the Z c states, suggesting a possible landscape of hadronic moleculear interpretations of the newly observed 
Z c states. We suggest that the 2 C (4200) meson is a DD* state while the 2 C (4020), 2 C (4050) and Z c ( 3900) mesons is 
either a D*D* or D\Di state. The stable DqD i molecular state does not occur in our result. 

The bottomonium-like qbbq molecular states are also studied and the numerical results are collected in Table El 
The extracted masses are predicted to be around 9.92 — 10.48 GeV, which are slightly lower than the masses of the 
charged 27(10610) meson. However, the charged 2f,(10610) meson is consistent with a BB* molecular state within 
the theoretical uncertainties. 

One finds that the hadron masses extracted from the color singlet-singlet currents are a bit higher than those 
extracted from the color octet-octet currents, for both the charmonium-like and bottomonium-like systems. This 
situation is different from the result in Ref. [49[, in which the color octet-octet tetraquarks were heavier. In Refs. [50, 
HU , the color-octet mechanism was found to give contributions to quarkonia production via the emission or absorption 
of a soft gluon in NRQCD. Similar mechanisms can be expected in the octet-octet quarkonium-like molecular systems, 
in which a color-octet Qq pair combines with another color-octet qQ pair by exchanging a gluon. 

The possible hadronic decay patterns of the qccq and qbbq molecular states can be discussed by considering the 
kinematic constraints and the conversations of parity, C-parity, isospin and G-parity. Considering the hadron masses 
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obtained in Table HI the possible S-wave two-meson hadronic decay channels for the charmonium-like qccq molecular 
states with I G J p ( G ) = are 


Z c DD*,D*D*,r lc (lS)p, J/4>n, 4>(2S)n , (23) 

and the possible P-wave decay channels are 

D*D, r) c (lS)h ((1235). (24) 

For their isoscalar qccq partners with I G J PC = 0 _ 1 H , the possible S-wave decay channels are 

Z c DD*,D*D*, Vc (1S)oj, J/tl’Vi J/H , (25) 

while the P-wave decay channels are 

Z c ->■ T7 c (lS’)fe 1 (1170), J/V>/o(980), J/^a 0 (980),/i c (lP)r ? . (26) 

For the bottomonium-like qbbq molecular states, the extracted masses in Table [TT] lie below the open-bottom thresh¬ 
olds so that only the hidden-flavor decay channels are kinematically allowed. The possible S-wave decay patterns for 
the qbbq states with I G J p ^ = l + l + (“) are Zb —> T(lS')7r, T(2S)tt while the S-wave decays are forbidden. For the 
isoscalar partners with I G J PG = 0”1 H , their possible S-wave and P-wave decay channels are Zb -A T(1 S)r] and 
Zb -A T(15)/ o (980), T(15)ao(980), respectively. 
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Appendix A: Expressions of spectral density for other interpolating currents 


In Eq. (TT71) . we have given the spectral density extracted from the current ■]{’. For other interpolating currents 
listed in Eq. ([5]) , we collect the expressions of the spectral density in this appendix up to dimension eight condensate, 
as shown in (fl5l) . 
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• For the current J. 
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P ( 2 )m (s) = 

P ( 2 1)(GG \s) = 
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• For the current J. 
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4a 3 /3 3 


TOq( 1 - a — /3) 2 (5 + a + /3) [(a + /3 )toq - a/?s]‘ 


„(8)<99) _ m Q(<7'7) 

P3 (Sj " _ ^ 4_ 

(9 2 s GG) 


(8)(GG), x 

P3 (S) = - 


9a 3 /3 3 

/■/tmax (1 — a — /3) [(a + /3 )toq — a/3s] [3wg(a + /3) — 7a/3s] 

da L. dp ~ -^- 

Wq(1 — a — /3) 2 (5 + a + /3) [mg(3a + 4/3) — 3a/3s] 


x m%n h'mtn 


x max / Hmax 


da 


dp- 


+ 


+ 


967t 6 J a . Jg . I 9a 3 /3 

" LX mtn h'mtn I ' 

7(1 — a — /3) 2 [(a + P)rriQ — aps ] [mg(13a + 13/3 + 29) + 12a/3s] 

288a 2 /3 2 

mg(l — a — /3) 2 [mg(o + /3) — 2a/3s] 
a 3 

7 [(a + /3 )toq — a/3s] [3mg(l + a + /3) + 4a/3s] 

48a/3 

(1 — a — /3) [(a + /3 )toq — a/3s] [mg(3 + a + /3) + 5a/3s] | 

12a 2 /3 J ’ 


4 8)< " G9> (8) = - 


mQ(gg s a ■ Gq) 
48?r 4 


da 


rfimax 

/ 

"Pmin I 


(1 — a — /3) [3mg(a + /3) — 5a/3s] 
a/3 


(1 — a — /3) [6wg(a + /3) — lla/3s] 71mg(a + /3) — 105 aps ^ 


2/3 




(8)(g q )(gG q ) () _ {qq)(qgs^-Gq) f 1 j 24m 4 Q (3 - a) 

P3 1 j “ 54m 2 J 0 2(1-a)a 2 1 q) 

TOg(24a 3 — 91a 2 + 54a + 9) 


2a(l - a) 2 


S (s - to|) + (13a + 12)37 (s - ™q) ( ■ 


For the current 


( 8 ) 


4m 


p[ 8)pert (s) = 


x max r h’max 


167T 6 


da 


x mtn " h'mtn 


d/3< 


(1 — a — /3)(1 + a + /3) [(a + /3)mg — a/3s] 
4a 3 /3 3 


TOg(l - a - /3) 2 (5 + a + /3) [(a + /3 )toq - a/3s] u 


pi 8)fe> W = 


2 tt 4 


(8) (GG) / X (ff s 2 GG) r— n 

Pi (s) =-/ da / 

" Olmin. " Pr 


96m 6 


9a 3 /3 3 

r/3 m ax (1 — a — /3) [(a + /3 )toq — a/3s] [3mg(a + /3) — 7a/3s] 

Ja 1 df> - ^ - 

TOg(l — a — /3) 2 (5 + a + /3) [mg(3a + 4/3) — 3a/3s] 
9a 3 /3 


Ot-min " Pmin 


dp- 


' CXmin " Pmin 


7(1 — a — /3) 2 [(a + P)rriQ — aps ] [mg(13a + 13/3 + 29) + 12a/3s] 

288a 2 /3 2 

mg(l — a — /3) 2 [mg (a + /3) — 2a/3s] 


(A6) 
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pf ){qGq \s) 


pf )m \s) 

p (8)(qq)(qGq >( fl ) 


7 [(a + /3)iriQ — a/3s] [3mg(l + a + /3) + 4a/3s] 

48a/3 

(1 — a — /3) [(a + /3)mg — a/3s] [mg(3 + a + /?) + 5a/3s] 1 


12a 2 /3 


m Q (qg s a ■ Gq ) 
48tt 4 


da 


rPmax 

/ 

“ ft min ^ 


7’ 

(1 — a — /3) [3m,g(a + /3) — 5a/3s] 


a/3 


(1 — a — /3) [6mg(a + ,0) — lla/3s] 71mg(a + 0) — 105a/3s 1 


2/3 


J’ 


10(2 


r nQ + s){qq } 2 


27t r 2 

(qq)(qg s cr ■ Gq) f 1 
54-7T 2 ,/n 


yJ\-4m 2 Q /s 


da< 


24m 4 (3 -a) . 

2(1 - a)a 2 5 


mg(24a 3 — 91a 2 + 54a + 9) 
2a(l - a) 2 


<5 (s — mg) + (13a + 12 )H (s — mg) > . 


(A7) 













